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Abstract — An  effective  Nonlinear  Schrodinger  Equation  for  propa¬ 
gation  is  derived  for  optical  dark  and  power  law  spatial  solitons  at 
the  subwavelength  with  a  surface  plasmonic  interaction.  Starting  with 
Maxwell’s  Nonlinear  Equations  a  model  is  proposed  for  TM  polarized 
type  spatial  solitons  on  a  metal  dielectric  interface.  Two  separate  sys¬ 
tems  are  considered  in  which  one  metal  dielectric  interface  has  a  dielec¬ 
tric  Kerr  medium  that  has  self-defocusing  and  another  similar  interface 
which  the  dielectric  Kerr  medium  that  has  self-focusing  depending  on 
the  modulus  of  the  electric  field  to  some  power  law  variable  p.  The 
beam  dynamics  are  analytically  studied  for  these  nano  waveguides. 

1.  INTRODUCTION 

Optical  spatial  solitons  have  been  studied  for  the  past  few  decades 
due  to  their  potential  applications  to  optical  pulse  compression  [4] ,  all 
optical  switching  [5],  logic  devices  [6]  and  others  [1].  The  advantage 
of  these  optical  beams  is  that  they  are  self  trapped  due  to  balance 
between  diffraction  and  nonlinearity  (self-focusing)  [16-19].  With 
this,  there  may  be  no  need  of  a  cladding  for  a  waveguide  to  keep 
the  light  beams  from  escaping  from  the  core.  These  spatial  solitons 
interact  with  one  another  while  displaying  properties  that  are  normally 
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associated  with  real  particles.  Bright  spatial  soliton  amplitude  profiles 
are  displayed  as  a  hyperbolic  secant  function  and  hold  there  profile 
at  short  distances  on  the  order  of  centimeters.  Dark  spatial  solitons 
are  similar  but  have  hyperbolic  tangent  type  profiles.  Dark  spatial 
solitons  were  first  observed  by  Jerominek  in  1985  and  Belanger  and 
Mathieu  in  1987  in  an  experiment  of  an  optical  branching  effect  in  a 
Titanium  and  Lithium  Niobate  planar  waveguide.  Spatial  solitons  have 
been  observed  also  in  other  experiments  [1-3] .  Gisin  and  Malomed  [7] 
showed  that  bright  and  dark  spatial  solitons  beam  size  cannot  be 
smaller  than  half  the  carrier  wavelength  in  the  transverse  magnetic 
(TM)  case.  In  the  transverse  electric  (TE)  case,  bright  and  dark  spatial 
solitons  are  subject  to  strong  instability  [8] .  Recently,  intense  study  has 
gone  into  an  area  of  nano-optics  called  plasmonics  due  to  its  potential 
application  for  development  of  smaller  optical  devices  such  as  nano¬ 
waveguides  circuits  at  the  subwavelength  [9—11].  Plasmons  were  first 
observed  as  quantized  bulk  oscillations  by  Pine  and  Bohrn  in  1952  as 
explanation  for  energy  losses  of  fast  electrons  passing  through  metal 
foils.  Surface  plasmons-polaritons  (SPPs)  are  a  quantized  excitation 
at  the  metal/dielectric  (or  metal/vacuum)  interface  with  one  material 
having  a  negative  permittivity  (usually  metal)  and  free  charge  carriers 
and  the  other  material  having  a  positive  permittivity  (dielectric  or  Kerr 
medium).  These  collective  oscillations  of  surface  charge  behave  like  a 
particle.  SPPs  support  phenomena  such  as  Surface  Enhanced  Raman 
Scattering  (SERS),  fluorescence,  enhanced  nanoscale  transmission,  and 
second  harmonic  generation  in  metals  from  subwavelength  slits.  These 
surface  plasmons  (SPs)  are  the  key  in  the  development  of  nano¬ 
antennas  and  platforms  for  sensing  biological  and  chemical  agents  [12— 
14].  Nonlinearity  plays  important  role  in  increasing  efficiency  of 
nanocircuit  and  waveguide  applications.  Past  studies  have  been  done 
on  nonlinear  guided  waves  in  metallic/dielectric  structures  localized 
in  the  transverse  dimension  [22,31].  Furthermore,  losses  in  the 
visible  frequency  range  have  to  be  considered  when  studying  soliton 
generation  and  propagation  for  plasmonic  structures  [1,2].  Davoyan  et 
al.  have  shown  that  the  bright  spatial  soliton  beam  size  limitation  can 
be  overcome  theoretically  by  using  surface  plasmons  [15] .  They  showed 
that  a  nonlinear  dielectric  depending  on  the  modulus  squared  of  the 
total  electric  field  and  metal  dielectric  interface  can  create  what  is 
called  “spatial  plasmon-polariton  solitons” .  This  was  verified  by  them 
through  full  vector  finite-difference  time  domain  (FDTD)  simulations 
with  commercial  software  Rsoft.  In  this  paper,  nonlinear  dielectrics 
and  the  behavior  of  SPs  with  self  defocusing  and  SPs  with  power  law 
nonlinearity  is  considered  theoretically  at  the  interface. 

In  Section  2,  a  nonlinear  Schrodinger  equation  is  derived  for 
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self-defocusing  based  on  metal-dielectric  boundary  conditions.  Next, 
the  one  soliton  solution  is  considered  from  the  effective  nonlinear 
Schrodinger  Equation.  Expressions  are  developed  for  the  integrals  of 
motion.  Spatial  frequency,  amplitude  relation  and  spatial  velocity  are 
also  derived.  Briefly,  the  requirements  for  dark  soliton  creation  and 
prediction  are  presented  which  apply  at  the  nanoscale  based  on  the 
effective  nonlinear  Schrodinger  Equation. 

In  Section  3,  an  effective  nonlinear  Schrodinger  equation  is 
considered  for  self-focusing  with  a  power  law  Kerr  like  medium  based 
on  metal-dielectric  boundary  conditions.  With  no  loss,  this  equation 
is  usually  solved  using  the  inverse  transform  method  [1,19].  In 
this  case,  a  secant  type  soliton  solution  is  proposed  and  conditions 
for  existence  are  derived.  The  results  include  equations  for  spatial 
frequency  with  a  nonlinear  correction  to  plasmon  frequency.  This 
also  includes  equations  for  the  amplitude  and  spatial  velocity.  In 
the  case  with  no  attenuation,  the  integrals  of  motion  are  conserved. 
Lastly,  the  case  with  attenuation  is  considered  which  the  integrals 
of  motion  energy,  momentum,  and  the  Hamiltonian  are  dynamic. 
The  integral  coefficients  in  the  proposed  model  nonlinear  evolution 
equation  contribute  to  derived  adiabatic  expressions  with  no  losses 
and  attenuation. 


2.  SELF  DEFOCUSING  WITH  METALLIC  DIELECTRIC 
INTERFACE 


In  order  to  create  self  defocusing  in  the  dielectric  nonlinear  medium  the 
linear  dielectric  constant  has  to  counter  the  intensity  of  the  dielectric 
medium.  Consider  an  illuminated  TM  polarized  light  that  is  created 
to  be  narrowly  focused  to  some  beamwidth  Rq  in  close  proximity  to  a 
microscopic  planar  waveguide  with  the  metal  on  top  and  the  nonlinear 
dielectric  on  the  bottom  (see  Figure  1  below).  Plasmons  are  excited  by 
light  scattered  on  the  metal  corner  so  plane  waves  will  not  be  excited 
in  the  dielectric  material.  In  addition,  a  thin  metal  screen  would 
have  to  be  place  in  front  of  the  dielectric  surface  [15].  To  derive  the 
ENLSE  (Effective  Nonlinear  Schrodinger  Equation)  we  start  Maxwell’s 
equations  for  planar  waveguide  interface  with  a  dielectric  and  a  metal 
(see  Figure  1  below).  But  first  we  consider  the  eigenmodes  for  the  case 
which  is  linear  (or  no  Kerr  effect  and  complex  dielectric).  Recall  that 
Maxwell’s  time  dependent  vector  equations  are: 


„  .  dB 

VxC+^“0’ 

„  ,  dD 

Vxh-w=0. 


B  =  noh, 

D  =  £0e  =  £0n2  C, 


(1) 
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Figure  1.  A  diagram  of  the  nonlinear  dielectric/nretal  interface  in 
order  to  avoid  plane  wave  stimulation  in  the  dielectric  material  [15]. 


where  yo,n,  and  £o  is  the  permeability  of  free  space,  the  index  of 
refraction,  and  the  permittivity  of  free  space.  Let  us  consider  £  and  h 
as  time  dependent  monochromatic  fields  defined  as 

h=  —  H(x,y,z)e~iut  and  C  =  E(x,  y,  z)e~iu}t.  (2) 

Hoc 

Substituting  Eq.  (2)  into  Eq.  (1)  reduces  to  time  independent  equations 

X7xE  —  ikH  =  0,  and  V  x  H+ikn2E  =  0,  with  k  =  2n/X.  (3) 

Taking  the  curl  of  the  above  Eq.  (3)  and  substituting  the  curl  of  H 
into  the  first  part  of  Eq.  (3)  and  normalize  by  1  jk  this  leads  to  the 
Wave  Equation: 

V  x  V  x  E  —  e(x)E  =  0.  (4) 

The  free  space  wavelength  is  A.  The  electric  field  E  is  the  new 
field  normalized  to  k.  In  this  case  assume  the  e  depends  on  x  for 
a  planar  structure.  Also,  we  assume  that  these  waves  are  transverse 
magnetic  (TM)  waves  so  E  =  (Ex,0,E~)  or  Ey  =  Hx  =  Hz  =  0.  Now, 
the  metal- dielectric  interface  is  designed  as  such  the  permittivity  for 
x  >  0,  e(x)  =  em  +  K  for  metal  and  for  x  <  0,  e(x)  =  Eun  ±  n\E\2  for  a 
Kerr  nonlinear  material  that  has  self-defocusing  (+  for  self- focusing) . 
S'  is  the  imaginary  part  of  the  dielectric  for  metal  and  r/  is  the  nonlinear 
coefficient.  Moreover,  the  dielectric  tensor  could  depend  on  either 
the  transverse  or  longitudinal  field  components  which  helps  simplify 
the  problem  and  has  been  verified  through  simulations  [32-34],  If  we 
separate  the  normalized  wave  equation  into  its  vector  components  the 


resulting  equations  are 

d2Ex  d2Ex 

d2E~ 

(5) 

dy 2  d z2 

for  x  and 

d2E-  d2E z 

d2Ex  ,  , 

for  z. 

(6) 

dx2  '  dy2 

-  dzdz*£<-x)E‘-° 
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Equations  (5)  and  (6)  cannot  be  solved  analytically  due  to  the 
dependences  of  field  components  to  nonlinear  propagation.  We  are 
looking  for  solutions  for  E  that  are  invariant  in  the  y  directions  so 

E  =  GExOe~i0z 

E-m^nl 
£m  T  £nl 

where  G  is  the  amplitude  and  j3  is  the  propagation  constant  or  plasmon 
wavenumber  normalize  to  k.  x  and  z  are  unit  vectors  in  the  x  and  z 
directions.  f5  is  also  the  effective  dielectric  constant  similar  to  effective 
index  of  refraction  usually  seen  in  classical  guided  mode  theory  [20]. 
Exo  and  Ezo  are  evanescent  type  fields  or  surface  plasmons  at  the  metal 
dielectric  interface  with  boundary  conditions 

Ex o  =  -^-eKlx  and  Ez o  =  eKlx  with  k\=02— em  for  x<0 

7  (9) 

Ex o  =— e  K2X  and  Ez o  =  e  K2X  with  k%  =  /32  — eni  for  x>0. 

Referring  to  the  above  Eq.  (7),  assume  that  TM  polarized  field  and  Ey 
is  neglected.  The  beam  in  the  y  direction  is  self-focusing  in  the  paraxial 
approximation.  These  are  considered  paraxial  plasmonic  beams.  We 
consider  a  solution, 

E  =  C\EXQe~i0zx  +  iC2EzOe~i0zz,  (10) 


x  +  iGEzOe~i0zz 
1/2 


(7) 

(8) 


with  C\  and  C2  being  slow  varying  amplitudes  and  Eq  =  (Ex o,  0,  Ez o) 
is  linear  Plasmon  profile.  Also,  the  plasmon  transverse  structure  does 
not  change  significantly  (or  dxC\  «  0  and  dxC2  «  0)  [15].  Substitute 
Eq.  (10)  into  Eqs.  (5)  and  (6)  taking  into  consideration  boundary 
conditions  (Eq.  (9)),  simplifies  to 


Exo 


d2c1  d2c1  dC! 

+  wrvr  -  2z/3 


dy2 


dz 2 


dz 


—  i- 


,dC2dEz  o 
dz  dx 


+kClEx0TV  \C!\2  E2x0  +  \C2\2  E2z0  C1Ex0  =  0,  (11) 


E 


z  0 


d2C2 
dy 2 


jdC'1  dEx.o 

dz  dx 


+kC2Ez0TV  \Ctf  E2x0  +  \C2\2  E2z0  C2Ez  o  =  0.  (12) 


Multiply  Eq.  (11)  by  Ex o  and  Eq.  (12)  Ezq.  Assume  C\  «  C2  with 
a  slow  vary  approximation  \d2C\/ dz2\  <C  \2i(3dC\/dz\  or  the  length 
scale  of  the  waveguide  (in  microns)  is  much  greater  than  the  carrier 
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wavelength.  With  neglecting  d2C\/ dz2 ,C\  =  C  and  integrating  with 
respect  to  x,  Eqs.  (11)  and  (12)  become 


—iP 


f  E20dx  dC  1  d2C  llvix) 


I 


Er 


+  -- 


1 

+*2 


fs(a 


dx 


Er 


dz  2  dy 


,2  T  2 " 


Er 


dx 


I 


Er 


\C\2C 


dx 


dx 


I 


Er 


-C  =  0. 


(13) 


dx 


A  more  familiar  form  of  Eq.  (13)  is 


with 


-l3W^  +  ]2W*^C  +  'nC=0 


w 

£ 

n 


f  ExOdx 


1 

2 

I 

E0 

dx 

f  V(x) 

Er 

4 

dx 

I 

fs(a 

EG 

0  - 

2 

i  t 

Eq 

ix 

2 

dx 

f 

Er 

2 

dx 

(14) 

(15) 

(16) 

(17) 


So  Eq.  (13)  is  a  nonlinear  Schrodinger  Equation  with  integral 
coefficients  that  are  a  function  of  x  for  a  dark  or  bright  solitons.  Again, 
rj(x)  and  $(x)  are  the  Kerr  nonlinear  term  and  the  imaginary  part  of 
the  dielectric  constant  for  the  metal.  Q  is  described  as  the  attenuation 
of  surface  plasmon  soliton  [15,  21,  22].  Considering  a  dark  soliton  and 
fl  =  0  or  no  attenuation  we  get 

-tl3W™  +  lw-m'2c=0-  (18) 

This  is  the  effective  nonlinear  Schrodinger  Equation  (ENLSE)  for 
plasmonic  dark  spatial  solitons  with  no  losses  with  an  integral 
coefficient  contribution  from  W  and  £  the  Kerr  coefficient  (x  <  0). 
With  £  =  0,  the  model  (Eq.  (18))  is  considered  a  model  of  beam 
diffraction  at  the  nanoscale  only. 

With  the  above  NLSE  equation,  there  exist  a  stationary 
fundamental  one  soliton  solution,  N  solitons,  and  even  algebraic 
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solutions.  For  simplicity,  we  consider  the  fundamental  mode  (or  lowest 
order  solitary  wave)  and  assume  a  stationary  solution  in  which  an 
Ansatz  is  described  as 

C(y,  z )  =  Ag[D(y  -  y)j  exp [i(-Ksy  +  usz  +  cr0)]  (19a) 

and  v  =  ^  (19b) 

dz 

which  the  function  g  is  described  as  the  shape  of  the  soliton  and.  u  is 
the  mean  velocity  of  the  soliton  with  mean  position  y(z).  The  constants 
A  and  D  represent  the  soliton  amplitude  and  width.  K,  los ,  and  <to  are 
the  soliton  frequency,  the  wave  number,  and  the  center  of  the  phase. 
This  frequency  uis  is  also  the  nonlinear  correction  to  the  plasmon  wave 
number.  The  phase  is  defined  as 

<f>(y,  z)  =  -ny  +  ujsz  +  <70 .  (20) 

With  r shape  =  D(y  —  y)  we  substitute  a  dark  spatial  soliton  solution 
for  g: 

C(y,  z)  =  Atanh(Tshape)  exp[i$(y,  z)\  (21) 

Substituting  Eq.  (21)  into  Eq.  (18)  we  get  the  following  relations 

k  =  (3Wv  (22a) 


Equation  (22e)  is  a  dark  spatial  plasmon  polariton  soliton  amplitude 
which  is  the  solution  to  the  ENLSE  Eq.  (18)  (see  Figure  2(b)). 
Dark  solitons  have  boundary  conditions  which  C(0,y)  approaches 
Co  the  background  amplitude  with  y  approaching  +oo  while  C(0,y) 
approaches  Coe10  as  y  approaching  — oo  with  theta  being  the  constant 
phase.  These  profiles  exhibit  an  intensity  dip  near  the  center  y  =  0. 
The  magnitude  of  the  dip  at  the  center  is  governed  by  the  amplitude 
A.  The  phase 


(23) 
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Power  Law  Soliton  (p=0.25)  Input  and  Output  Beam  Profile  and  Spectrum  with  p=1  initial  profile 


Normalized  space  t 


Normalized  Frequency  to/2*n 

(b) 


Figure  2.  This  is  Split  Step  Fourier  simulation  of  the  normalized 
NLSE  Eq.  (14)  with  =  0  for  a  soliton  calculating  power  and  spectral 
power  normalized  amplitude  ijj  with  initial  amplitudes,  (a)  Power  law 
bright  soliton  p  =  0.25  with  p  =  1  initial  profile  for  comparison,  (b) 
0.8  for  a  dark  soliton  without  attenuation. 

is  a  single  parameter  that  characterizes  the  dark  soliton  solu¬ 
tion  [1,  3, 19].  The  dark  spatial  plasmon  polariton  soliton  can  be  cre¬ 
ated  with  an  arbitrary  small  dip  on  a  continuous  wave  background. 
The  condition  for  existence  of  a  dark  soliton  at  the  nanoscale  would 
be  that  the  input  beam  would  be  [1,  3] 

C(0,y)  =  C0  exp(i6>)  +  C01(y)  (with  Cm  0  as  |y|  ->•  oo)  (24) 
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and 


A  =  Re  ^exp(-i0)  J  Cm(y)dy  j  <  0. 


(25) 


Under  this  condition  there  exist  two  discrete  eigenvalues  in  the 
associated  inverse  scattering  problem  such  that 

7i,2  =  ±Coi(l  -  1/2A2)  (26) 

This  corresponds  to  a  two  dark  solitons  with  equal  amplitudes  and 
opposite  velocities  [1,3].  With  the  above  conditions,  the  intensity  dip 
always  produces  a  pair  of  dark  solitons. 

The  equations  for  the  integrals  of  motion  renormalized  due  to 
concern  of  divergence  are 


Pr  =  /  (Cq  -  \C\2)dy 


Mr=  ~  (CyC  —  CyC*)  (  1  - 


— oo 

oo 


C l 

|C|2 


dy 


(27a) 

(27b) 


Hr  = 


V 


dC 


dy 


|C|2 


+  /  [F{Cl)-F{I)}dl\dy  (27c) 


c20 


with  F(I)  =  —I  and  F(Cq)  =  —  Cq  for  the  dark  plasmon  soliton. 
Pri  Mr,  and  Hr  are  renormalized  soliton  power,  momentum,  and  the 
Hamiltonian  quantities.  So  the  equations  for  power  and  Hamiltonian 
for  the  dark  plasmon  soliton  are 

Pr  =  2  Vi  A  (28) 

Hr  =  4  ^(C2-V2)l.  (29) 

To  include  attenuation,  Eq.  (14)  for  self-defocusing  is  normalized  to 
form  the  effective  NLSE 

\*r*+**-°  <3»> 

with 


C  Ro 


z  —Rq 


— /31U\  1/3  , 

—7—  1  v  y  = 

-/3W\1/3 


Ro  6 


e  ) 


■a  p  = 


(-PW)  '2 
QRl 

2  (-/31U)2/3 


e  ) 

(0_1/3 


(31) 


92 


Crutcher  and  Osei 


amplitude  =  075  Input  and  Output  Beam  Shape  and  Spectrum  1-secht2 


(a) 


(b) 


(c) 

Figure  3.  With  initial  profile  0)  =  1  —  -Bosech2(r)  simulation 
of  the  normalized  NLSE  overcome  by  losses  with  calculation  of  (a) 
power,  (b)  spectral  power  and  (c)  contour  plot  of  the  beam  for  a  dark 
plasmon  soliton  with  attenuation  em  =  —25  +  il .74,  £unear  =  4.84, 
Rq  =  0.35  pm,  rj  =  3.0  x  10~6  V/m2  and  Bq  =  0.65. 


where  Rq  is  the  initial  beam  width.  Again, with  P  =  0,  Eq.  (30)  is  the 
normalized  effective  NLSE  which  has  solution  of  hyperbolic  tangent 
function  form  (See  Figure  3).  Integrating  W,  14,  and  £  which  are 
indefinite  and  assuming  rj  and  ?  are  constants,  the  coefficients  depend 
on  dielectric  constants  for  the  metal  and  Kerr  material.  There  is 
a  discontinuity  at  the  origin  as  far  as  the  boundary  conditions  (see 
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Figure  4.  (a)  Initial  profile  -0(7”,  0)  =  1  —  -Bosech2(r)  simulation 

plasmon  soliton  with  Bq  =  0.2  without  attenuation.  An  unstable  dark 
spatial  plasmon  soliton  simulation  showing  break  up  into  Gaussian 
like  spatial  profiles,  (b)  0)  =  -Botanh(r)  with  Bo  =  1.  (c) 

V>(r,  0)  =  1  —  l?osech2(r)  with  Bq  =  0.8. 


Eq.  (9)).  In  applying  the  Split  Step  Fourier  Method  Figures  3  and  4 
below  describes  the  propagation  if  P  /  0  with  the  initial  profile  a 
hyperbolic  tangent  with  em  =  — 25+*1.74,  £unear  =  4.48,  Ro  =  0.35  pm, 
rj  =  3.0  x  10“6V/m2,  and  amplitude  0.8.  The  dark  plasmon  soliton 
breaks  up  into  small  Gaussian  like  spatial  profiles  moving  into  the  ±r 
directions.  With  tj  and  ?  depending  on  x,  the  integral  coefficients 
profiles  would  have  to  be  determined  for  dielectric  and  nonlinear 
materials  for  certain  frequencies  which  are  not  considered  in  this  paper 
due  complexity  of  the  integral  coefficients. 

It  has  been  shown  that  dark  spatial  solitons  exist  when  attenuation 
is  ignored.  The  soliton  features  at  the  nanoscale  are  similar  to  optical 
spatial  solitons  above  the  one  half  carrier  wavelength  limit  with  the 
appropriate  assumptions  and  approximations.  The  attenuation  has  to 
be  managed  in  order  achieve  self  defocusing  and  generate  dark  solitons. 
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3.  POWER  LAW  MEDIUM  DIELECTRIC  METALLIC 
INTERFACE 


Spatial  solitons  have  been  known  to  be  created  in  materials  with 
power  nonlinearity  such  as  semiconductors  and  other  materials  [23,  24]. 
Considering  dielectric  s(x)  is  not  dependent  on  the  square  of  the 
modulus  of  electric  field  but  the  electric  field  to  some  number  say  p, 
the  equation  for  the  dielectric  is: 

S(x)  =  Sunear  +  T] \E\2'P  (32) 


The  power  law  in  NLSEs  has  been  studied  extensively  for  optical 
temporal  solitons  [18,19,25-28].  For  p  =  1,  Davoyan  et  al., 
theoretically  and  numerically  using  FDTD  (Finite  Difference  Time 
Domain) ,  showed  bright  spatial  plasmon  polariton  soliton  propagation 
for  a  nanosized  dielectric  interface  [15].  Following  the  same 
assumptions  for  Eqs.  (1)-(10)  similar  in  the  previous  section,  but 
use  Eq.  (32)  instead  of  a  Kerr  nonlinear  material  (p  =  1)  based  on 
Maxwell’s  equations  we  get: 


Ex0 

+  T] 

Ezo 

+  T] 


^d2C1  d2C1  dCx 
+  wrw  - 

dz 


2 p  JT^P 


dy 2  dz2 
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(33) 


dz 


ICi] 


2 p  z 7>2p 


Ko  +  ICjI 


dx 

2 P  p]2P 


zO 


C2Ezq  =  0, 


Again,  we  assume  Eq.  (33)  has  the  same  boundary  conditions 
as  Eq.  (11),  with  E0  =  (Ex0, 0,  Ez0),  E  =  CiEx0e~tf3zx  + 

iC2Ezoe~lPzz,  dxC\  «  0  and  dxC2  «  0,  Ci  ~  C2  =  C,  and  the  slow 
varying  envelope  approximation,  \d2C\/ dz2\  <C  \2i(3dC\/ dz\  is  applied. 
Multiplying  the  first  equation  by  Ex o  and  the  second  by  Ez o  taking  the 
sum  and  integrating  over  transverse  component  x,  we  get  the  power 
law  effective  nonlinear  Schrodinger  equation: 
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BC  1  B2C 

-i(5W—  +  -  — -v  +  &p\C\2pC  +  inC  =  0  with  0  <  p  <  2.  (35) 
Bz  2  By 1 

iid 

W_  fExOdx  ,aRn\ 


f  Eq  dx 


f  Tf(x)  E%  +  E%  [E20  +  E2z0]  dx 


— *  " 

f  Eq  dx 


l  f  s(x)  E0  dx 

n  =  - - - .  (36c) 

f  Eq  dx 

p  is  limited  to  0  <  p  <  2  due  to  self  focusing  singularity.  Taking  into 
consideration  a  generalized  NLSE  such  as 


l2pVv  =  o, 


.Bl/Jg  1  (  32^g  32^g  \  ,  ,  2p  ,  _n  ,  Q, R, ^ 

*' +  2  V +w  (36d) 

with  Dim  being  the  transverse  dimension,  there  are  3  different  cases 
to  consider:  1)  In  the  critical  case  p  ■  Dim .  =  2  solutions  become 
singular  in  a  finite  2.  2)  With  p  ■  Dim ,  <  2,  which  is  the  subcritical 
case,  the  diffraction  always  dominates.  3)  In  the  supercritical  case 
p  ■  Dim  >  2,  there  is  a  large  class  of  smooth  initial  amplitudes  for 
which  a  focusing  singularity  forms  in  a  finite  distance  z.  With  the  case 
between  the  subcritical  and  supercritical  self-focusing,  near  balance  is 
achieved  between  focusing  nonlinearity  and  diffraction.  With  =  0, 
Eq.  (35)  the  model  is  described  as  diffraction  at  the  nanoscale.  We  can 
solve  for  a  stationary  soliton  solution  by  making  II  =  0  for  attenuation. 
We  propose  the  following  Ansatz  [29,  30]: 

C{V’ 2)  =  coshHBfo  -  jj)]  eXp(-“!' +  +  ia)  (37) 

with  v  =  g. 

Here,  B,  A,u,K,  and  us  are  the  width,  amplitude  of  the  soliton, 
mean  velocity,  spatial  frequency  and  the  wave  number,  a  and  y  is 
centre  of  the  soliton  phase  and  the  centre  of  the  soliton.  Substituting 
Eq.  (37)  into  Eq.  (35)  with  Q  =  0  using  the  relation  m(2p  + 1)  =  m  + 1 
for  powers  of  hyperbolic  cosine  function  in  the  resulting  relations  are 

k  =  (3Wv  (38a) 


2p2f3W 


(38b) 
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and 


1 
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C(y,z) 


A 

T  “ 

coshr  [B(y  —  vz—  y)\ 


exp  (—iKy+iuz+ia). 


(38c) 

(38d) 


(e) 

Figure  5.  Simulation  of  power  law  plasmon  soliton  with  i/j  = 
i?osech1/,p(r),  Bq  =  0.8,  p  =  0.25.  (a)  With  no  amplitude  dependence 
in  the  phase,  (b)  With  small  attenuation  ^  0.  (c)  Phase  depend  on 
amplitude  with  fi  ^  0.  (d)  Stable  with  proper  amplitude  dependent 
phase  (Eq.  (38c)).  (e)  0  /  0  but  complex. 
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The  case  p  =  1  reduces  to  the  Kerr  law  for  the  bright  spatial 
plasmon  soliton  [15].  In  this  case  our  soliton  parameter  dynamics  are: 


dB  „  d.A  dy  k  dn 

—  =  0  —  =  0  —  =  -  — 

dz  dz  dz  /3W  dz 


(39) 


The  integral  coefficients  affect  the  power  law  soliton  propagation, 
but  do  not  change  the  usual  properties  (see  Figure  5).  The  three 
conserved  quantities  energy  (E),  momentum  (M),  and  the  Hamiltonian 
(H)  for  the  power  law  are 
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with 
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where  the  Gamma  function  is  defined  as 
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(43) 

(44a) 
. .  (44b) 


The  conservative  quantities  depend  on  the  initial  beam  profile  and 
the  integral  coefficients.  Referring  back  to  Eq.  (35),  with  14  7^  0  the 
conservative  quantities  would  change  along  with  dynamic  variables.  To 
slowly  add  loss  to  the  Eq.  (35)  we  multiply  a  <  1  to  the  loss  term. 
The  new  equation  with  a  perturbation  is 

Pt( 7  1 

-if5W—  +  -^+^P\C\2pC  +  iaQC  =  0  (45) 

If  we  normalized  Eq.  (45),  then  resulting  equation  is 


‘aS'+12w+^2Pc'--‘£c‘ 

with  transformations  defined  as 
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With  the  loss  term  include  the  adiabatic  variation  of  the  soliton 
parameters  are 
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Transforming  back  from  B  to  z  from  the  original  NLSE  and  calculating 
integral  on  the  right  hand  side  of  Eq.  (47)  through  (49) 
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Equation  (50)  through  (52)  is  ordinary  (and  possibly  coupled) 
differential  equations.  The  change  in  the  soliton  k  and  amplitudes 
ASzB  depend  on  conserved  energy,  attenuation,  and  the  integral 
coefficients  of  the  initial  input  beam.  Also,  the  change  in  energy, 
momentum,  phase,  and  velocity  as  the  evolution  of  the  center  of  mass 
are: 


j  jy  _ 

—  =(0W)-12s/^rnaE 

4=0 

dz 


(53) 

(54) 


“  =  p=&r(pwy1 


B2  Io, 0,2,0  K 

1 

j 

2  70,2,0,0  2 

7q,2,0,0 

A2p 


2  2 

coshr  (s) 


ds 


(55) 


100 


Crutcher  and  Osei 


V 


dy 

dz 


bp(PW)-1 


a 

Io,2, 0,0 


2  Ban^&p  f 

A2  J 


(56) 


00  2 

with  J7>5)VlW=  /  v^g<(v)(^y(^)™dv  and  u  =  B(y  -y(z)). 

—  OO 

The  function  g{v)  is  the  profile  of  the  ansatz  solution.  With 
the  loss  term  included  in  equation,  one  could  also  obtain  a  quasi¬ 
stationary  solution  but  this  will  not  be  included  here.  Here,  with  some 
assumptions  it  has  been  shown  that  it  is  theoretically  possible  to  have 
power  law  spatial  soliton  propagation  at  the  nanoscale  due  to  surface 
plasmons. 


4.  CONCLUSIONS 

Although  spatial  solitons  are  diffraction  limited  in  a  Kerr  medium 
alone,  this  limitation  can  be  overcome  by  a  nonlinear  dielec¬ 
tric/metallic  interface.  It  has  been  shown  for  this  interface  with  a 
self  defocusing  medium  dark  spatial  solitons  exist  theoretically  called 
dark  spatial  plasmon  polariton  solitons.  The  proposed  model  for  these 
dark  plasmon  solitons  was  a  nonlinear  evolution  equation  called  the 
effective  (1  +  1)  dimensional  (first  1  is  beam  spread  and  the  second 
1  is  propagation  coordinate)  Nonlinear  Schrodinger  Equation.  Using 
soliton  theory,  adiabatic  parameters  were  derived  with  the  attenuation. 
Another  type  of  interface  was  considered  applying  a  power  law  material 
with  a  metal  and  then  a  power  law  type  (1  +  1)  dimensional  Nonlinear 
Schrodinger  Equation  was  derived.  Lastly,  the  attenuation  was  consid¬ 
ered  in  the  power  law  derived  adiabatic  equations  for  changes  in  the 
energy  or  power,  momentum,  spatial  frequency,  amplitude,  velocity 
and  phase. 

One  could  also  derive  Algebraic  and  N  soliton  solutions  for  these 
dark  and  power  law  spatial  plasmon  solitons.  One  may  also  consider 
an  approach  with  dielectric/metallic  planar  waveguides  existence  of 
vector  solitons  or  derivation  of  a  (2  +  1)  dimensional  NLSE  model.  A 
more  accurate  approach  to  the  above  two  systems  would  be  to  consider 
numerically  the  coupled  Nonlinear  Schrodinger  Equation  like  equations 
which  removes  some  approximations  and  assumptions.  The  dark  and 
power  law  spatial  plasmon  solitons  attenuate  over  short  distance  as 
they  propagate  along  the  dielectric  metallic  interface.  This  can  possibly 
be  minimized  by  considering  solitons  in  tapered  plasmonic  waveguides, 
a  subject  of  future  investigation. 
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